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An analysis^ Is developed that maKCs it possible to 
integrate von Karnian^s boundary- layer momentum equation 
directly. For this purpose the skln-fr let ion coeffi- 
cient is expressed as a function of the Reynolds nuinber 
based on the boundary-layer momentum thickness and the 
boundary -layer shape is aas-uuned to bo constant. The 
inte.e;rsted equation permits the bo;.ndary- layer moment-am 
thickness at the airfoil trailing edge to be rapidly 
deteririned for the estimation of airfoil profile-drag 
coef f c3 ent s . 



IrriRODUCTIOIT 



'"ecent increases in size and speed of modern air- 
craft have emphasized the need for a method of esti- 
mating rapidly the profile drag of wings. One of the 
most tim.e-ccns"Uining operations in a profile-drag estimate 
has been the calculation- of the thickness of the turbu- 
lent boundary layer. Squire and Young (reference 1) 
calculate the thicknesc of the t^irbulent boundary layer 
by a step-by-step integration of von E'arman's boundary- 
layer irom_entum equation (reference 2). A logarithmic 
skin-friction formula designed to fit Schlichting ' s ap- 
proximate formula for the skin friction (reference 5) is 
used in the integration. Reference li gives a method for 
the determination of the burbulent boiindary- layer thick- 
ness vvhich eliminates the step-by-step computations of 
reference 1 but requires that suitable approximations to 
the airfoil pressure distribution be made. The method of 
reference S eliminates both the step-by-s"Lep computations 
of reference 1 and the necessity for making approximations 
to the airfoil pressure distribution required in reference 
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however, reference 5 restricted to a particular 
boundary- layer shrpe and does not explicitly indicate 
how to use an arbitrary sjcin-friction law with the 
inte.r^rated von Karmah equation. 

In the present paper the step-by-step coirputation of 
reference 1 is eliminated along with the necessity for 
riaklnr approximations to the airfoil pressure distribu- 
tions, which were required in reference Ij.. A more general 
specification of the boundary-layer sliape than that given 
in reference 5 i.? used and a method for use of any skin- 
friction lav' with the integrated von Karman e^quation is 
presented. The Integrated form of the von RaVmp'n momentum 
equation appears applicable either to the laminar or the 
turbulent boundary la^ror provided that the appropriate and 
constant shape of the boundary layer i.'^' assumed and the 
appropriate formula for the skin friction is used. For 
the tiu^bulent boundary layer, a skin-friction formula 
based directly on the data obtained from experiments in 
pipes (reference 2) .is presented in a form readiljr usable 
for boundary- layer caD.culations • 

'Hhe present work provides a method for the calcula- 
tion of the thickness of the turbulent boundary layer 
which is moi'e rapid than that of reference 1, more flexible 
than that of reference and more general than that of 
reference 5* The flexibility of the method arises from 
the fact that the effects on the calculated turbulent 
boundary- layer thickne^.v^ of local changes in the airfoil 
pres3iu'-e distribution, changes in the boundary- layer shape, 
and changes in the skin-fr ir^t ion lorm.ulas can all be 
easily eva]_uated. The method, however, like those of 
references 1, i|, and 5 cannot predict the onset of turbu- 
lent separation and, if the danger of turbulent separation 
is thought to exist, the method of reference 6 may be used 
to check for such a condition. Although the present method 
may be used in such cases to obtain miOm.ent^am thickness up 
to the point of separation, it should jiot b e used in the 
dete-^mination of the drag coefficient si The present mxethod 
can be used to calculate momentiom thicknesses for flow 
over rough as well as smooth surfaces, if a curve of skin- 
friction coefficient for the rough surface is available. 

The skin-friction equations presented should provide 
a good estimate of the turbulent boundary- layer thickness 
for use in profile -drag estimates when the airfoil is 
smooth and there Is no danger of turbulent separation. 
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SYIvIBCLS 

c airfoil chord 

Cr\ S3cticn urofile-drag coefficient 

K boundary-layer shape paraineter f or v/hich values 

are given in reference 6 {j^J 

k constant in equation for sl-in-fr ict ion coefficient 

n con^.tant in e-iju&tion for skin-fr lotion coefficient 

q dynamic pressure, pounds per square foot (^f^} 

R^. wing Reynolds number based on chord and free- 

stream velocity 

Rq boundary- layer Reynolds number [ij""^ 

s distance along airfoil suz'face, usually taken 

equal to x 

s^ position on surface at start of integration 

u local velocity inside boundary layer 

IT lo^.al velocity outSide boundary layer 

free-stream velocity 
TJ-j^ velocity st Sq/^- 

X distance alonr airfoil chord 

y perpendicular distanc: from wing siu'face 

5 full boundary-layer thickness 

5"' displacement thickness \ f (l-|j)dy 

9^ momentum thickness at Sq/c 

skin friction per unit area, pounds per square foot 

u kinem.atic viscosity 

The subscript t is used to indicate the condition at the 
trailing; edgie . 
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p density 
9 miomentum thickness 
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ATTALYSI3 



The von Karman boiindary-layor monentuun equation, 
which is used to cslculrte the Tnor'ientum thickness of the 
boundary layer, is given in reference 2 and may be vjr it- 
ten as 

ds 2 \qlds 2q 



The mor:.entui7i equation is a differential equation that 
can oe integrated if the bovindary-layer shape, and there- 
fore H, is assumed to be constant and if the local skin- 



friction coefficient 
function of R9. 



2q 



is expressible as a power 



In order to integrate the equation the skin-friction 
coefficient is written as 



k _ kU' 



n 



where k and n are constants that are to be determined. 
The mO'Ttentum equation is vTitten as 



de ^ H 4- 2 dU g 

ds U ds 



ku^ Q-n 



rjn 



v/hich is a differential equation of the Bernoulli type. 
After the equation is into^:;:rat6d and the arbitrary con- 
stant arising from the integration is evaluated, the 
result is 



IT 



Oy 



.11+2 



s / c 
n+1 , 



(l+nik l-^^^ I IT 



/cV^o 



IT- 



l+n)k / 
p n J, 

(II+2) (n+1) 




(K + 1 ) (n+l)+l 



1+n 



(1) 
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The definite integral occiirrlng In equation (1) may be 
evaluated by eibher an analytical or a graphical metliod, 
v/hichever is more convenient o 



DISCUS3I0JI 



Vtoen the rnomentim equation w£?s integrated, the 
assuiniotion was made that the shape of the boundary layer, 
and therefore H, was constant over the chord. Al- 
though the assumption is not true in general, the error 
in 9 will usually be small because of the manner in 
which H appears in the moment-^jmi equation. 

H]xperience has indicated that, for the calculation 
of the thickness of the tm-^bulent boundary layer, the 
value of H may be chosen close to l.L for airfoils 
that have sm.all pressui^e recoveries but, for airfoils 
that have high pressure recoveries, v/hich are thick or 
operating at high lift coefficients, H should be chosen 
closer to 1.6. In refereiices 1, Ij., and 5, H was chosen 
as constant and equal to l.i.,i. The thickness of the 
laminar bo^undary layer can usually be obtained with suf- 
ficient accuT'acy for profile-drag calculations if the 
integrated momentum equation is used, together with a 
value of H of 2.^92, which is the value for the Blasius 
fist-plate profile (reference 2). 

In order to use the present form of the integrated 
momentum equation, it is necessary that the skin-friction 
coefficient be written as 



2q o n 
^ R9 



For the lam.inar boundary layer having the Blasius shape, 
n = 1 and k = 0.2205. 

A skin-friction formula for the tui-^bulent boundary 
layer which has the required form is that of Falkner 
(reference 7) v/ho, after analyzing various data obtained 
from flat plates, suggested the relation 

ISL = 0>006555 

2q R 1/i 
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A skin-friction forr:ula based directly on pipe experi- 
ments can be derived oj use of equations (r!!2.2l), 
(22. lo), and (22.17) reference 2, pa^es I42 and 1)4.5. 

This forn.ula 



is -^Dlotted in figiu-^e 1. The cuj-^ve of equation (2) and 
the'relatlon suggested by Falkner (reference 7) are 
plotted in figure 2. At values of Rq above about 2000, 
the ci-ir "'&s are fo->md to be in good a£;reerrient but at 
values below 2000 diverge considerably. 

The skin-friction formula derived in reference 3 
indicates skin-friction coefficients' at small values of 
Rq that are in good agreeinent with equation (2). The 
skin-friction coefficients given by Falkner «s equation 
are tlxerefore probably too small at small values of Rg . 

Equation (2) is based directly on tk-e pipe experi- 
ments to v.hich the formula of Squire and vcung (refer- 
ence.. 1), although not derived in a direct manner, also 
ov/es its origin. This equation and a convenient plot 
of the equation are given in reference 60 Consideration 
of flat-plate skin-friction data and approximations, 
which were made in the interest of simplicity of expres- 
sion, account for the difference between the Squire and 
Youjag formula and equation (2). The skin-friction coef- 
ficients indicated by use of the Squire and Young formula 
are slightly lov;er than those obtained by use of equa- 
tion (2) at large values of Rg but are about the same 
as those obtained by equation (2) at small values of Rg. 
The introduction of equation (2) also serves to indicate 
how an arbitrary sV:in-fr ic t ion formula miay be used with 
equation (1) . 

In order to permit the integration of the mom:entum 
equation, the skin-frict ion relation chosen for use 
is approxim.ated by a sim.ple power fujiction of Rg . This 
approximation is m_ade by using two values of Rg and the 
skin-.friction coefficients associated with the two Rg 
values to evaluate the two constants, n and k, in the 
equation 





1 



2 



(^) 




1^ 



n 
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n 



The 1 omnia s to ''oe usod are 

R9-|_ 



k 



where the subscript 1 identifies the skin-friction 
coefficient occurring at Rq^, the smaller value of R9, 

and the subscript 2 denotes the skin-friction coefficient 
occurring at Pg^, the larger value of Rg. Experience 

has indicated that, if Rg^ is chosen about ten times as 
larre as , th^; aoor oximation to ociuation (2) will be 

good over the entire range of Rg values of the turbu- 
lent bomdary layer likely to be encoujitered in a spe- 
cific case. The beginning of the turbulent boundary 
layer has been foLind to be the region in v/hich it is most 
important to have a good approximation of the skin- 
friction ciH-'ve. Because of this fact, Rg.^^ should be 

Chosen close to the initial value of Rg of the turbu- 
lent boundary layer, vrhich is ass-umed equal to the final 
value of Rg^ in the laminar boundary layer. 

After the momentum thickness of the boundary layer 
at the trailing edge has been obtained, the profile-drag 
coefficient is calculated by the relation developed by 
squire and Young in reference 1 

H+5 




2^ 



The subscript t is used to identify the condition at 
the trailing edge. 
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COMPARISON WITH EXPERII'iEITT 



Calculations of profile-drag coefficients have been 
made by the method deGcribod herein and the results have 
been coinpared with those obtained by other investigators. 
The profile-drag coefficient of the 25-percent-thick 
airfoil at a lift coefficient of O.25 and a Reynolds 
nvtmber of £.2 x 10°, for which data are given in refer- 
ence 1, was calculated and compared with the experimental 
value of 0.0080 given in reference 1. iTnen the equation 



v/as used, the calculated pro ile-drac coefficient 
was 0.0076. "'lien the equation 



v/as used to fit equation (2) for values of R9 between 
1,100 and 11,000, the calculated drag coefficient 
was 0.0077. Squire and Young in reference 1 obtained 
a calculated value of O.OO79. In all three cases the 
value of H v/as taken equal to l.i;. 

calculations of the boundary- layer momentum thick- 
nesses, along the chord and the profile-drag coefficients 
of the NACA '0012 airfoil at zero angle of attack for 
four Reynolds niombers have been made and compared with 
the experimental results obtained from reference 8. 
The comparison is shown in the following table: 



"^o _ 0.0063.33 



2q ' Re^/^ 



1° - 0'0090^ 

2q ~ Rq 0.2028 




Experi- 'Calculated, Calculated, 
mental | Squire and Falkner's 
(ref- lYoujig (ref- equation 



Calculated 
equation 
(2) 



2.675 X 10 i 0.0071 

5.780 1 .0070 

5.550 .0068 

7.560 ' .0067 




0 . 007ii. 
.0072 
.0071 
.0071 



0.0067 
.0069 
.0069 
.0069 



0.0069 
.0070 
.0070 
.0069 
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The profile-drag coeff Iclents that were calculated by 
using Falkner's formula in one calculation and the 

relation 



I2. - ^*0093k 

2q Pq0.206S 

as an approximation to equation (2) "betv/een Rg values 
of 900 ^"^'^ 9000 in the second calculation are ^iven in 
this table^ together with the experimental values from 
reference 8 and profile-drag coefficients calculated by 
the squire and Young method in reference 8. The 
boimdary-layer momentun: thicknesses over the surface of 
the M7iCA 0012 airfoil ^ which were calculated by using 
the approximation to equation (2) and the formula sug- 
gested by Falkner, are given in figure 3 together with 
the exf)erlmental values of iromentijun thiclcness obtained 
from reference 8. The differences between the two cal- 
culated momentum thicknesses are caused by the use of 
the two different skin-friction fcriiiulas, the value of 
H being taken as I.5 for both calculations. 

In order to indicate the effect of the choice of H 
for the turbulent boundary layer on the calculated 
profiler-drag coefficient of a thin airfoil at a small 
lift coefficient, the profile-drag coefficient of the 
NAOA 0012 airfoil at zero angle of attack and a Reynolds 

number of T.SbO x 10^ was also calculated for values 
of H of I.Il and 1.6. The profile-drag coefficients 
were O.OO68 for a value of H of l.k and O.OO7O for a 
value of H of 1.6. The change in H caused a rela- 
tively small change in the profile-drag coefficient of 
the NACA 0012 airfoil at zero angle of attack because 
for this airfoil at zero angle of attack the change in 
velocity over the airfoil surface is small. As can be 
seen from, equation (1), a change in the value of H 
will have greater effect on the calculated profile-drag 
coefficient for cases in which large changes in velocity 
occm'' (for example, on thick airfoils or on airfoils at 
large lift coefficients) than for cases in vv^hich the 
change in velocity is small. 
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SAMPLE CALCTJLATIO:i 



A calculation is given for the profile-drag coef- 
ficient of the MACA 0012 airfoil at a Reynolds nuiiroer 

of 2.675 X 10^ and zero angle of attack. The velocity 
distribution, v/hich was obtained from reference 8, is 
r;iven in figure ii and is the experimental velocity dis- 
trioution, uncorrected for tunnel-wall effects. At the 

Reynolds number of 2.675 ^ 10^ the transition point 

was at — = O.I4.8. 



The equation that was used for the calculation of 
the morrent-umi thickness of the laminar boundary layer is 
obtained from equation (1) by placing; ^ 0/0 - 0 , n-1, 
k = 0.2205, and Fi = 2.592. The resulting equation is 



c 





8.18]+ 




q.66Il 



d^ 
c 



Vs/c 





1 

2 



Then 




c 



The final value of Rq in the laninar boundary, which 
is assurried to be the initial value of Rq of the 
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tiirbulent bo-undary layer, is then calculated: 



f=o.i,8 



= 2.675 >^ 10^ X O.COO29OI X l.ll^i. 

- 888 



A value of 9OO was chosen for Rg^ and 9OOO for Rg^ 

for the appi'o.xlrnation of eq.uation (2) by a power fuaction 
of Rg . The computations for n and k are as fol- 
io v,-s : I^ror,^ figure (1) 



,^ = QOO 



~ = 0.002289 
2q 



= 9000 



T 

— = 0.001421 



n = 



1r = 



.. 0. 002.2 89 
'^^o.ooiJ.i2i _ o.j}.76 

log 10 2.502 

0.2063 

0.002289(900)'-^*^'-''^^ 



Tnen, 



= O.OO93I1. 



I2. = O.Q093li- 
2q p 0.2068 
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The forrnula for the turoulent portion of the 
boundary layer then h6coi.ies , for H - 1.5^ 



c 



r 0.01126 /jj V4- • 017^ s 



■t- ( 0 . 00029 01 ) ^ • ^ ( 1 . ll;ii ) ' ^ • ''-"'^ 



1/1.207 



The valii.c of 

'0 



I (Jl.] d- Is eaual to O.676 



8.i8ii 




T"i fifrui'e ^ are shov/n ciirves o.f [ — j and 

'ttJi..017 ' ^ VV 

— ] that an-oly to the lar-'lnar and turbulent 

Po/ 

portions of the boundary layer, respectively. Then 
= 0,002005 for one s'oi'face and 

^■1 

C^^ J: X 0.002005 X 0.5545-25 z: 0.0069 

CONCLUSIONS 

1. If the skin-friotion coefficient is expressed as 
a pov;er function of the Reynolds number based on the 
boundary- layer momenturi thickness and the boundary- layer 
shape IS assumed to be consta:it, it is possible to inte- 
grate von Karman's boundary-layer nomentum equation 
directly;- . 

2o The integrated boundary- layer momentum equation 
permits the boundary- layer momentum thic'icness at the 
airfoil trailing edge to be ranidly determined for esti- 
mation of airfoil profile-drag coefficients. 
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Figure 2.- Plot of Palkner»s equation and equation (2) against Rg . 
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(a) R, 2.675 X 10 . 
Figure 5»- Experimental and calculated momentum thickness on 
the NACA 0012 airfoil. 
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(b) R, 3.780 X 10^. 
Figure 3.- Continued. 
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(c) R, 5.350 X 10^. 
Figure 5.- Continued. 
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(d) R, 7.560 X 106. 
Pig\B*e 3*~ Concluded. 
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Pigxire i^..- Velocity distribution for NACA 0012 airfoil at 

(Taken from reference 8.) 
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Fig. 5 
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\ (H + l)(n + 1) +1 
Figure 5.- Cvtrves of (■^j for the 

\ 0/ 

NilCA 0012 airfoil at a = 0°. 
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